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In this paper we study and compare susy unification using two different ap- 
proaches in order to take into account the effect of fight particle thresholds on the 
evolution of gauge couplings: the step-function approximation, on the one hand, and 
a mass dependent procedure, which gives a more accurate description of the depen- 
dence of the results on the masses, on the other. We also include the effect of heavy 
thresholds, when SU(5) is chosen as the unifying group. We find that the mass- 
dependent procedure excludes scenarios where all susy masses are below 1 TeV, and 
favors a value of a^{mz) near its upper experimental bound, contrary to the results 
obtained with the step-function approximation. We underline the dependence of the 
results on the procedure chosen to deal with light thresholds. 
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Experimental data on sin 2 9 and a e measured at rriz can be used to test supersymmetric 
grand unification theories (SGUT). Use of 2-loop RGE's to go from Mj (the unification 
scale) down to mz, gives a value for sin 2 9 in good agreement with experiment []3J. It can also 
be checked that when the three gauge couplings ,ai2 and a^ 1 of SU(3) x SU(2) x U(l) 
run from low to high energies, they become equal at a scale Mx high enough to satisfy 
experimental lower limits on proton decay f|]||. Further work has been done on this problem, 
making use not only of limits on proton decay [§],[|] , but of cosmological arguments on the 
relic abundance of the light susy particle (LSP) || to constrain the susy spectrum. In 
general, the susy spectrum obtained from this scenario is below 1 TeV, and available in 
the next generation of particle accelerators. Their detection would be the best test of susy 
theories. 

In studying SGUT's one has to deal with both light thresholds (associated with the susy 
masses), and heavy thresholds (associated with the heavy masses) arising from the specific 
unification group G. When supersymmetry is broken at the Planck scale by a "hidden" 
sector, the large number of susy masses can be determined at the weak scale in terms of a 
small number of parameters at the unification scale: mi/ 2 (universal gaugino masses), m 
(universal scalar masses), A (cubic scalar couplings) and m 4 (bilinear scalar coefficient) [|7|]. 
On the other hand, for the minimal choice G = SU(5), there are three basic heavy mass 
parameters: My (heavy gauge bosons), M$ (heavy colored scalars), and M s (heavy adjoint 
scalar multiplet). At issue is how to handle these thresholds. 

In general, thresholds are included in the evolution of the gauge couplings by a step- 
function, given logarithmic corrections. For light thresholds, instead of using specific susy 
masses one can used an effective scale M susy to summarize the effect of the degenerate spec- 
trum ||. Corrections due to heavy thresholds in general increase M susy , so that taking a 
null correction will give a lower bound on M susy . Nevertheless, heavy threshold corrections 
are constrained if one takes into account limits on proton decay via dimension-five opera- 
tors (allowed in SGUT's, via the $ field) ||, and the theoretical requirement that Yukawa 
couplings of heavy scalars do not blow up below the Planck scale [f| . 
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In this paper we use a different approach (different from the step-function approxima- 
tion) to include threshold effects in the study of SGUT's. In a recent paper || we have 
defined effective charges [[UJ for the three gauge couplings, calculated with a mass depen- 
dent subtraction procedure , which include a complete dependence on the masses. With 
this method, one includes, for example, threshold contributions due to massive gauge bosons 
W ± and Z°, which are missed in the step-function approximation. In the step-function ap- 
proximation each particle with mass m > mz contributes to the running coupling a^ 1 ^) 
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with a logarithmic term In on the other hand, the mass dependent method gives a more 
precise description of the physics via a function, for both m > mz and m < mz, which can 
be approximated by 



fin, m ) = ln —f-. r ' W 

H l + c m 2 



where "c" is a constant of order 1-10. If m 2 <C m\ we recover the usual term ln^J-, and 
when m 2 ^> fi 2 there are no contributions: the particle decouples from the theory [12 



The constant "c" is chosen so that it "matches" the logarithm with the exact function for 
intermediate scales, mz < m < /i, and its value depends on the type of particle running 
inside the loop. When one has two mass degenerated fermions or two scalars, "c" can be 
obtained from the leading term of a power expansion of the exact threshold function when 
m 2 //i 2 goes to infinity. In other cases, the coefficient given by the power expansion has to 
be slightly modified to get a best fitting for intermediate scales. 

To see the basic features of the mass dependent procedure versus the step-function, we 
first study the evolution of the coupling constants from mz to Mx without any reference to 
heavy thresholds. Thus, we impose the unification condition: 

aj\M x ) = a^(M x ) = a^(M x ) = a~ G l , (2) 

with Mx > 10 16 GeV to prevent fast proton decay. The effective charges a^ x (Mx), 
a^Mj), aJ~ 1 (Mx) calculated at one-loop order are given by: 

ar\M x ) = a-\mz) + ±- £ bf ] f k \M x , m k ) , (3) 
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where we sum over all the particles in the SSM, and f {k) {M x , rrik) is given by ([!]). The values 
of rriw and m z are well known, and we have as arbitrary mass parameters: m t (top quark)[], 
and nih (light Higgs particle of the SM); (winos) and m g (gluinos); rriq (squarks) and mr, 
rrif (sleptons left and right); m M (higgsinos); mu , itl+ and m a (scalar masses from the second 
Higgs doublet needed by susy). For the susy masses we take a simplified parametrization 
in terms of and £o = (^o/ m i/2) 2 ? neglecting the mixing between winos and higgsinos, 
and the s-top left and right: 

m a = 7774/2 , m- g = ?>m 1 /2 



nip, 



m i/2\f7 + £0 , (for all squarks, including s-top) 



m 



L 



777-1/2^0. 5 + £ , 777^= 7771/2^0.15 + ^0 



we also take 777 + ~ 7?7 a ~ mjj. Thus, we are left with six arbitrary parameters: m t , m^, 



m i/2i £o> m ii and 777//. These parameters have lower bounds [14] derived from experimental 
searches for the top, Higgs and susy particles: 



m t > 91 GeV , m h > 60 GeV |L3 
7771/2 > 45 GeV . 

Moreover, perturbative bounds on Yukawa top couplings and quartic Higgs couplings yield 
the theoretical upper bound m t , rrih < 200 GeV ||15|| : and no extreme fine-tunning on the 
susy parameters gives 7771/2, tt7o < 1 TeV. 

We also need the initial values OL[ x {m z ). The values for af^m^) and a2 1 (m z ) are 
derived from the experimental data on sin 2 # and a" 1 , 

a-\m z ) = 127.9 ±0.2 |§, sin 2 6(m z ) = 0.2327 ± 0.0007 [|TJ, 



but for a^{m z ) there is no agreement between different measurements [I8|: The latest LEP 



data average to a^{m z ) = 0.122, while data from low energy measurements average to 



The remaining fermions in the Standard Model are taken as massless. 
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Oiz{mz) = 0.109. Because of this discrepancy we will not take a^{mz) as the initial data, 
instead we will derive it from the unification condition. In this way susy masses can be 
bounded by requiring that a^mz) be in the range (0.108, 0.125). 

In Fig. (1) we have plotted a^m?) versus log 10 (mi/ 2 ), for different values of the remain- 
ing free parameters, and including the experimental and theoretical constraints on m^, a% 1 
and Mx mentioned above: 

TOi /2 > 45 GeV , M x > 10 16 GeV , 0.108 < a 3 (m z ) < 0.125 . (4) 

Since m M , itlh contribute with the same sign, we have simply taken m M = itlh in order to 
check the basic features of the mass dependent method. Furthermore, since both aj 1 (m^) 
and M x depend very slightly on £ , for this plot we have fixed £ — 1- Varying the susy 
masses, we observe a trend similar to what happens when using a step-function: the higher 
the susy masses, the higher a^ l {mz) and the lower Mx- However, while with the step- 
function (Fig. 2) one obtains the limit a 3 (m z ) < 0.116 (this procedure does not distinguish 
masses lower than mz), with the mass dependent procedure the full range of a 3 _1 (m^) can be 
covered for adequate values of the mass parameters. Thus the limiting values a^(m z ) = 0.125 
and Mx = 10 16 give us a lower bound on m M = m# and an upper limit on m^. These 
bounds depend on m t , and £o; both of them decrease with m t and m^, while the bound 
on mi/2 increases with £ and the bound on decreases, remaining practically unchanged 
f° r Co > 10 4 . Therefore, if we allow a maximum difference of two orders of magnitude 
between m and mi/ 2 , we get the following upper bounds (using the central values of a e and 
sin^) 

m 1/2 <2.5TeV (m h = 60GeV, m t = 9lGeV), 

> 338 GeV (m h = m t = 200 GeV , m^ = m H ). (5) 

The one-loop calculation already shows the differences between the step-function and 
the mass dependent procedure. The first favors a^imz) to be in the range of the low energy 
experimental data; in addition, the lower data excludes susy masses greater 1 TeV. On the 
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other hand, when we include a more precise treatment of thresholds, the naturalness bound 
of 1 TeV favors a^{mz) to be in the range of the last LEP data. 

When we improve the accuracy, calculating at two-loop order, the value of 013 (mz) 
diminishes by about 10%. In the case of the step-function, this effect drives the upper 
bound obtained at one-loop towards the experimental upper value, while the naturalness 
bound on susy masses gives a lower bound a^{mz) > 0.1180. With the mass dependent 
procedure, a 10% decrease in the one-loop result puts the values of a 3 1 (rriz) obtained with 
= ran = 1 TeV away from the experimental band, being now the lower bound on 

= rriH around 10 TeV. In Fig. (3) we have plotted the upper bound on m.1/2 and the 
lower bound on = rtin obtained at two-loop order. The behavior of these bounds differs 
slightly from one-loop: the upper bound on ra\M is a maximum for £ — 60 (the value 
plotted), nearly independently of m t and m^, and the lower bound on m M decreases even 
for £0 > 10 4 . We would need £0 > 10 10 to have m M in the range of TeV, and then squark 
and slepton masses of order of PeV . Thus, this simple scenario of perturbative unification 
at two-loop order is not compatible with experimental data on a 3 (m^) and the theoretical 
naturalness bound on the susy spectrum. We need at least a heavy higgsino or a heavy 
Higgs (or both of them) beyond this bound. We also note that the remaining susy masses 
(gauginos, squarks and sleptons) are allowed to have values below 1 TeV . 

Up to now we have not included the effects of heavy thresholds, but a correct picture of 
perturbative unification needs to include them. Thus, the unification condition @ reverts 
to: 



2 As was pointed in Ref. [19|, this bound is consistent with the one obtained including the impact 
of the evolution of gaugino masses (EGM) in the step-function approximation pQ] ; the lower bound 
on a^{mz) (upper bound on (mz)) comes from taking mw 2 ~ 45GeV, for which the EGM 
effect is small. This does not occur with the "mass dependent" effect, which gives differences of 
order 13%, no matter the region of masses we take. 
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47r«: 1 (fi) = 4vraG 1 (/i) + A;(/x, Mj) , (6) 

where Mj are the heavy masses, and /x is a mass scale satisfying <C /x <C Mj, i.e., 
it is far away from both light and heavy thresholds [pl] , p2| . The unification condition is 
obtained taking into account the decoupling of the heavy degrees of freedom from the low 
energy theory; integrating out these fields from the action one gets an effective field theory 
in terms of low energy parameters (c^, mj), which are related to the high energy parameters 
(«g, Mj) through Eq. (f|). The functions Aj at one-loop include logarithmic terms due to 



heavy degrees of freedom, and constant terms due to light degrees of freedom, |2l 



A,,( / x,M,)=Af ) + A^ ) (/i,M J ) , (7) 

A (//, Mj = 7cj In -T H In -T F In . 9 

fi 3 fi 3 fi 

(s = 2.029884...) 

Here h, f refer to light scalars and fermions respectively, and H, F to the heavy ones; c$ = 
C 2 (Gi) for d = SU(3), SU(2), 17(1), and cj = C 2 (G) - C 2 (Gi); and T* are representation- 
dependent co efficient s0. Thus, if we make the reasonable assumption that the members of 
each heavy supermultiplet are degenerate we get: 

n/r on \ i 

(10) 



Al(yU, Mj) 


66 

= y + 


96 , 
— in 
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My 4 
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(11) 
(12) 

The couplings a~ l (fi) in Eq.(|) are the same as given by (3), where Mx is now replaced 
by fi <C Mx- Since heavy masses are typically of order 10 16 GeV , and light masses are 



3 The constant term — Cj/3 is not present in when one uses the DR subtraction procedure 
[23], as it occurs in supersymmetric theories [24]. 
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expected to be less than ITeV, we choose fj, = 10 7 GeV. Eliminating a G 1 (//) from (^), we 
obtain a^ 1 (mz) and In My: 

1 13 3 

ou^imz) = - (3aJ 1 (mz) - a^K)) + —^3(7714, ^) - — lnM$ - — , (13) 

2 v ' 87T 07T 07T 
37T 3 

In My = — \ Zai X {m z ) - a^ l {m z )) + — Fy(m;,/i) 

_ A i n M „ - - In M E - — + . (14) 
40 * 8 L 10 8^ K J 

The dependence of a 3 _1 (m2) and In My with the light masses (given by the functions 

F 3 (nii,fi), Fy(rrii,}i)) is qualitatively the same as before, i.e., (mz) increases with the 

susy mass parameters while In My decreases. We now focus our attention on the heavy mass 

parameters M$ and M E . We see from Eq. ( |i~3"D that (mz) depends only on M$, so the 

limits on M$ will put a bound on 0:3 (m^). The lower bound comes from the experimental 

limits on proton decay via dimension-five operators. Using a chiral Lagrangian technique, 

the lifetime obtained for the dominant mode is 0: 



r(p -> K+Vp) = 6.9 x 10 



31 



0.003 sin 2p H M # 10~ 3 2 



yrs , (15) 



13 1 10 17 /K~,m r , 77^) 
where yet three more unknown parameters have popped-in: the hadron matrix element 
parameter (3, which ranges from 0.003 to 0.03 GeV; the ratio of vacuum expectation values 
of two Higgs doublets tan/3#; and the parameter y tK , which represents the ratio of the 
contribution of the third generation relative to the second. To allow an M$ as low as 
possible, we take (3 = 0.003, sin2/?# = 1, ancf] 1 + y tK = 1. The experimental limit for 
this mode is r(p — > K + u^) > 1.0 x 10 32 yrs [p5] , so we get: 

M$ > 1.2 x 10 20 f(m q , m h mm) = M% in . (16) 

The function f(rriq, mj, m^,), with the parametrization we have adopted for the susy masses, 
is given by: 



In this case we take the simplest choice, because of our experimental ignorance about the value 
of y tK ■ In fact, y tK could be negative, giving |1 + y tK \ < 1. 
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/K, m h me) = — ^— ( ln(£o + 7) - ln(£ + 0.5)) , (17) 

6.5mi /2 \ £o + 6 6-0.5 / 

and therefore the lower bound only depends on mi/2 and £o, decreasing with m\/2 and m . 

On the other hand, both M$ and Ms can be bounded from above by requiring that the 
Yukawa couplings involving these fields do not blow up below the Planck scale ||. This 
leads to M$ < 2M V , M s < 1.8M V , and from (|TJP we get the upper bound on M$ in terms 
of the light masses as well as Ms: 

157T / i i \ 15 . 

lnM$ < — (a^ [m z ) - a 2 (m z ) ) + —F v (mi,n) 

5 , 52 200s 40, , _ mnT , . 

lnM s + = + — h2 = hMT • 18 

37 37 111V3 37 * v ; 

With M™ n (m!/ 2 , 6), M$ ax (mi, M s ) we get upper and lower bounds on a^" 1 (m^), which 
have to be within the range of experimental data. Therefore, we can play with the expres- 
sions fll3|), (|i~6|) ) (|T8|), and the limits on the susy masses trying to check whether or not the 
perturbative unification scenario is compatible with all the constraints. 

In order to compare, we first examine the case of the step-function approximation for 
F 3 (mj, /i) and Fy(rrii, y). With the constraint of susy masses below 1 TeV there is no problem 
in having a 3 inside its experimental range. As it can be seen in Fig. (4), the upper limit 
on £q gives us the minimum allowed value for M$ and the maximum for M^; furthermore, 
the lower limit on Ms would give us the lower allowed value for £ , which is just reached 
for the maximum allowed values of and m , and the upper one for M$. In principle, 
we do not have any constraint on the lower bound for Ms, except the requirement that 
there is no large splitting between the heavy masses. As M$, My will be around 10 16 GeV, 
we use Ms > 10 13 GeV to give the results. Increasing Ms reduces both M^ ax and a™ ax . 
These results are summarized in Table I, where we have considered limiting susy masses of 
ITeV and 2 TeV respectively. For both examples, the allowed range on a^{mz) is inside 
the experimental bounds; raising the limiting mass enlarges it, the same happens with 
the allowed range for M$ and m^. It is seen that if we require all light masses below 
ITeV, and Ms > 10 13 GeV, we get mus < mz, near it lower experimental bound, and also 

9 



strong constraints for m , M$, M s , so that practically, m ~ ITeV, M$ « 10 16 - 6 GeV, and 
M s w 10 13 GeV. There are no further constraints on m^, m H except, of course, the common 
upper limit we chose for the susy masses. 

As we have printed out without taking into account heavy thresholds, the mass dependent 
procedure gives us values of a^ l {m z ) lower than those obtained with the step-function 
approximation. And now, to satisfy the constraints on M$ and a^irnz) we will need to 
have mo, or ran beyond ITeV. The lower £ (^o)j the higher will be M$, and higher 
values of m M , ran will be needed to raise a^ l (mz) above its minimum bound (Fig. 5). The 
results obtained in this case are given in Table II: there are no solutions if we maintain all 
susy masses below 1 TeV . 

We have seen that the use of a mass dependent procedure and the requirement of hav- 
ing susy masses not too high, favor a value of a^ l (m z ) (a^[m z )) near its lower (upper) 
experimental bound. This is easily understood, since the mass dependent procedure has 
the effect of raising the values of a^ l (^) at high energies, with respect to those obtained 
with other methods (MS, step-function), and their different results tend to merge as susy 
masses are increased. When the maximum allowed value for a^ l (mz) is not enough to 
unify the couplings we will have to increase the masses; sometimes above its naturalness 
bound, as it happens in the two-loop calculation without heavy thresholds. Improving our 
scenario of coupling constant unification with the effect of heavy thresholds does not im- 
prove the situation about the light masses. In principle, taking M$ low enough we will get 
ctz(mz) < 0.125. But experimental data on proton decay do not allow to freely decrease 
the value of M$ without adequately increasing the squark and slepton masses. Thus, in the 
end, we conclude that it would be necessary to have some of the susy light masses, (either 
squarks and sleptons, or higgsinos or heavy Higgs), heavier than the naturalness bound of 
1 TeV commonly taken, contrary to the results obtained using the simplest treatment of 
the light thresholds given by the step-function approximation. 

We do not want to emphasize the numerical results, which depend on details of the model 
(such as susy mass parametrization) as well as on the experimental data. We would specially 
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like to advert to their dependence on the procedure chosen to compute the thresholds, which 
becomes relevant when running the couplings over such a huge range of scales: from m z 
to Mx (all light thresholds are crossed). It is clear that the step-function approximation 
gets worse as the number of crossed thresholds begins to proliferate; this is the case, for 
example, in susy theories, where the mass dependent procedure gives the exact contribution 
for each massive degree of freedom, independently of their number, and avoids the presence 
of cumulative errors that afflicts threshold crossing in a less complete treatment. Therefore, 
in models with intermediate mass scales or with many heavy-matter degrees of freedom, the 
use of mass dependent procedures is mandatory before reaching conclusions on unification 
and allied phenomena. 
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FIGURES 

FIG. 1. Values of a^ 1 (mz), compatible with the unification condition (Eq. [2]), calculated 
with a mass dependent procedure at 1-loop order, for different values of = ran: rnz, 1 TeV, 
WTeV, 100 TeV. Dotted lines are the experimental limits on oi^ (rnz) and mj/ 2 ; solid lines are 
for mt = rrih = 200 GeV, and dashed lines for rat = 91 GeV and = 60 GeV . The straight lines 
(solid for ra t = = 200 GeV and dashed for rat = 91 GeV, rrih = 60 GeV) are the upper limit 
obtained for a^ 1 (mz) when imposing Mx = 10 16 GeV (m^ = ran increase along these lines from 
bottom to top). The allowed region for a^ 1 (mz) are to the left of the straight lines, and between 
the dotted lines, 8 < a^ 1 (mz) < 9.2 and m 1 / 2 > 45 GeV. 

FIG. 2. Same as Fig. [j], but with a^ 1 (mz) calculated with the step-function approximation. 

FIG. 3. Lower bound on = ran (bottom curves), and upper bound on rnwj (straight 
lines), calculated with a mass dependent procedure at 2-loop order. The lower bound on is 
obtained for the lower experimental limit on ql^ (rnz) and the upper limit on £o (f° r the plot we 
choose £o = 10 4 , see text). The upper bound on m\i 2 is obtained for the upper limit on Mx, and 
£o — 60. Solid and dashed lines follow the same convention as Fig. [j] and ^. 

FIG. 4. a^ 1 (mz) versus log 10 (M$) at 1-loop order with the step-function approximation, for 
limiting susy masses of 1 TeV, ra 1 / 2 = 45 GeV (solid lines) and m™^ = 53 GeV (dashed lines), and 
satisfying different constraints: (i) The less slope bottom lines are obtained with M™ n (mi/2j£o) 
(Eq. JO]), and fixing = ran = mz] £o increases along these lines, from right to left. We have 
marked the points $?™(mffi = 45 GeV), and ^(mf/f = 53 GeV). Dotted lines set Mf in 
and M$ ax . (ii) The most slope lines are obtained with M£ Lax (m 1 / 2 , £o, m M ) (Eq. |18|) , with £ fixed 
(really depend very slightly on this parameter), = ran increasing from bottom to top, and 
limiting values for M s , Mg in = 10 13 GeV to the right, and Mg ax = W 13A GeV to the left. 
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FIG. 5. Same as Fig. ^j, now with a mass dependent procedure and limiting susy masses of 
2 TeV (there is no solution for 1 TeV); m™^™ = 45 GeV in solid lines and = 86 GeV in dashed 

lines. Dotted lines for M™ n , M^ ax , and lower experimental value for a^ 1 (mz) (a^ 1 (mz) > 8) . 
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TABLES 



^max 

m 


tmax 
SO 


M min 


M™ ax 


a 3 


m l/2 


tmin 
SO 


m min 


M min 




a 3 


1000 


494 


1Q16.58 


1Q 13.4 


0.112 


53 


357 


936 


10 13 


1Q16.63 


0.119 


2000 


1975 


1Q16.06 


10 16.2 


0.108 


221 


82 


936 


10 13 


igl6.63 


0.122 



TABLE I. 1-Loop with Step-function" 

a Mass values in GeV 



^max 

m 


tmax 
SO 


m min 


M min 


M max 


_ min 
a 3 


^max 
m l/2 


train 
SO 


M min 


M max 


a 3 


1492 


1100 


1492 


igl6.29 


10 16.34 


0.125 


45 


1100 


10 13 


igl6.29 


0.125 


2000 


1975 


805 


igl6.06 


10 16.37 


0.123 


86 


540 


10 13 


igl6.39 


0.125 


3000 


4444 


358 


10 15.75 


igl6.40 


0.120 


333 


81 


10 13 


igl6.54 


0.125 


4000 


7901 


179 


-j_q15.53 


10 16.41 


0.118 


4000 





10 13 


igl6.64 


0.125 


5000 


12348 


106 


-j_q15.35 


10 16.42 


0.116 


5000 





10 13 


10 16.72 


0.125 


6084 


18278 


m z 


10 15.20 


1Q 16.44 


115 


6084 





10 13 


10 16.79 


0.125 



TABLE II. 1-Loop with MDSP' 



a Mass values in GeV 
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